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Abstract 

Presently, the main methods for describing a non-equilibrium charge-transporting steady state 
are based on time-evolving it from the initial zero-current situation. An alternative class of theories 
would give the statistical non-equilibrium density operator from principles of statistical mechanics, 
in a spirit close to Gibbs ensembles for equilibrium systems, leading to a variational principle for the 
non-equilibrium steady state. We discuss the existing attempts to achieve this using the maximum 
entropy principle based on constraining the average current. We show that the current-constrained 
theories result in a zero induced drop in electrostatic potential, so that such ensembles cannot 
correspond to the time-evolved density matrix, unless left- and right-going scattering states are 
mutually incoherent 
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I. INTRODUCTION 



The problem of charge transport through nanoscale objects became a very intensive area 
of research during the last few years. Numerical calculations reflect the rapidly advancing 
development of experimental techniques exploring typical characteristics of atomic wires 
or single-molecule junctions^. The most usual theoretical methods used to address these 
issues are either direct occupation of scattering states^ or using the non-equilibrium Green's 
function formalism for steady-state (NEGF)^. Both of these rely upon an assumption 
that the use of a ground-state density-functional theory gives the effective self-consistent 
potential. This assumption is now known not to be correct even though corrections to it 
might be relatively small, which is particularly true for systems with open channels^. The 
NEGF method makes use of semi-infinite reservoirs which demand further approximations 
to the effective potential if one is to be able to perform ab initio studies of the systems of 
interest, e.g. typically the electronic structure of the leads are taken as that of the leads in 
equilibrium. The errors included this way can be assessed and calculations for representative 
systems suggest them to be no more than a few percent^. A second source of discrepancies 
between the experimental and theoretical work comes from the ambiguity of the geometry of 
the molecular junction one employs. The most common choice is the ground state geometry. 
Even though several interesting results concerning non-equilibrium forces exis t 10 ! 11 , ab initio 
studies of current-distorted geometries are still absent. 

Almost all of the above-mentioned development is based on time-evolution as a means 
to obtain the non-equilibrium steady stat o 12 i 13 i 14 , although in most cases the time evolution 
is used only formally for the derivation of the formulae used within the NEGF formalism. 
The other alternative would be to use some sort of variational principle that would directly 
lead us to the non-equilibrium density matrix (or state). While it might be difficult to be- 
lieve that such a principle would exist for a general dissipative system, its existence for the 
purposes of quantum transport is relatively easy to accept. At the level of a non- interacting 
or mean-field description of the electrons, it has been known for some time that such a 
variational principle exist a 15 i 16 , and it can be extended formally to fully interacting elec- 
tron systems, provided that there is a physical relaxation process that "washes out" initial 
correlations^. These are based on the search for a state with minimum energy, consistent 
with prescribed numbers of right- and of left-going electrons (hereafter referred to as the 



2 



"left-right" scheme, LRS) or prescribed average current (current constrained approaches). 
The LRS prescription is directly motivated by the Landauer-Biittiker formalism^ developed 
for mesoscopic systems and can be generalized by maximizing the information entropy with 
constraints on the average energy (corresponding to the introduction of the temperature 
in equilibrium) and average numbers of left- and right-going particles To remove the con- 
cept of single-particle orbitals from this approach, Frensley suggested to constrain the local 
Wigner distribution functions for electrons with positive and negative momenta in the left 
and right lead respectively^. This has the advantage that it can be used also for interacting 
approaches as has been demonstrated recently by Delaney and Greer^. A different exten- 
sion of the single-particle theory was presented by HershfielcU^, who formally constructed a 
non-equilibrium steady-state density matrix using many-body field operators corresponding 
to generalizations of the scattering states. 

As a alternative to the density constraint, several authors have suggested the total cur- 
rent as a means to keep the steady state out of equilibriu m 15 i 21 i 22 i 23 i 24 i 25 i 26 i 27 . The current 
constraint, unlike the density constraint, has immediately a well-defined form even for in- 
teracting electrons. It was soon realized that the current-constrained density matrix corre- 
sponds to a situation different to that obtained within the LRS, although no clear consensus 
exists regarding what experimental situation it describes. It has been mentioned already in 
the work of Ng^ that the theory should correspond to a const ant- current experiment (as 
opposed to constant- voltage). However, this cannot be the only criterion, since the system 
of interest is very small and so the ability of the system to explore the whole Hilbert space of 
admissible density matrices should be considered with care. The latter property is indeed at 
the center of the formulation employing constrained searches within the maximum entropy 
principl e 26 ' 28 . In this respect it is interesting to note recent work by Di Ventra and Todorov 2 ^ 
who have formally considered a variational principle based on constraining the current for 
a quasi-steady state of discharge of a large but finite electrodes through a nanoj unction. 
The steady-state-current constrained ensembles, which are the subject of the present paper, 
should correspond to a long-time and infinite-size limit of their considerations. 

The various current-constraint formulations differ in some details. Ngi^ considers a treat- 
ment where the current operator is altered in such a way that it does not mix the right- 
and left- going scattering states. To achieve this one has to drop all the off-diagonal matrix 
elements of the current operator in the scattering-states representation. This eventually 



3 



leads to a theory that is similar to the LRS, but with the occupancies now depending on 
the current that the particular state carries, as well as its energy. The applied bias is deter- 
mined from the difference in the local electrochemical potential between the left and right 
asymptotic regions which, strictly speaking, corresponds to the electrostatic drop around the 
sample. Heinonen and Johnso n 21 ! 22 consider only systems that are translationally invariant 
along the current flow. Under such circumstances the off-diagonal matrix elements of the 
current operator are absent by symmetry and the notion of applied bias is only formal. They 
determine the latter from the analogy with the scattering-states-occupation theory and the 
resulting I — V characteristics in the linear regime are identical to the LRS description. One 
should point out that this analogy fails as soon as one enters a strongly non-linear I — V 
regime with, for example, a current flow through a resonant barrier. 

In both of these treatments the effective Hamiltonian with the current constraint com- 
mutes with the physical Hamiltonian or, in other words, that the constrained density matrix 
is stationary. The approach developed by Koso v 23 ' 24 departs from this point and instead 
constrains the current to be uniform throughout the system. As a result, the density ma- 
trix is not time-independent which brings into question its relevance for the description of 
a steady state. The problem of the steady-state character of the density matrix obtained 
from a current-constrained search has been studied by the authors of this paper—. It has 
been found that the steady-state requirement does remove most of the off-diagonal matrix 
elements of the current operator, as is assumed by Ng. However, those off-diagonal elements 
between states with the same energy do not disappear, which leads to a density matrix 
different from that anticipated by Ng. The induced drop in the potential was found from 
the local neutrality conditions in the asymptotic regions, in a manner similar to that of Ng, 
and the paper additionally discussed the relation of the Lagrange multiplier A (that imposes 
the current constraint) to the applied bias voltage. 

The aim of this paper is to present a clear relationship between the above-mentioned 
current-constraining schemes, as well as to discuss their limitations. We start by analyzing 
the role of the external applied bias within the formalism. This yields a link between the 
time-evolved and variational approaches, and gives a supportive argument for the form of 
the steady-state requirements implemented in our earlier work^. In Section IHII we address 
the induced drop in the electrostatic potential. Conversely to what has been claimed before 
we have found that unless we remove all of the off-diagonal elements of the current operator, 
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the induced drop is exactly zero. This holds for the uniform-current approach, discussed in 
the Section where all elements are retained 2 ^ as well as for the approach (used in our 
previous study^) in which the current operator has off-diagonals only between equal-energy 
states. In Section II VI we describe how many-body interactions of the electrons with the 
environment correct the previous results to a physically meaningful picture where nonzero 
drop in the potential is found and the agreement with the conventional approaches (NEGF, 
LRS) achieved within the linear response regime. 



II. APPLIED EXTERNAL FIELD AND MAXIMUM ENTROPY 



We describe the steady-state-current situation using the maximum entropy principle for 
the statistical density operator p2£ 



5 \ -Tr [p log(p)] - TrA * [Mil \ = 0, (1) 



where A, are Lagrange multipliers that guarantee chosen values, Ai, of the averages of chosen 
operators A*. That is, instead of following the specific time evolution caused by an applied 
external field, we assume that the final steady-state is representable by the statistical density 
operator p that possess the same current and/or induced drop in potential. While all of the 
existing schemes agree on using the total energy and total number of electrons as two of the 
constrained operators Aj, the constraint that keeps the system out of equilibrium varies and 
is either the total current or the occupancies of right- and left-going scattering states. 

In real-time evolution approaches, there are no non-equilibrium constraints or multipliers. 
Instead, the applied external field acts as a driving force for the current flow. The character 
of the external field seems to be a source of certain confusion in the community. Some authors 
use a ramp-like external potential that has a finite slope between contact o 20 ' 24 even within 
maximum entropy schemes, while some avoid its presence altogethe r 15 i 17 i 21 i 22 i 26 . On the other 
hand, in calculations based on the time evolution (scattering states, NEGF formalism), the 
role of the applied external field is frequently circumvented by the application of a difference 
in electrochemical potential between two initially isolated leads. The electrons are taken 
to be non-interacting while in the leads. This construction, however, leads to violation of 
local charge neutrality in the non-interacting regions as discussed elsewhere^. Neither does 
it clarify the relation between the total and induced electric field. Clearly, the problem with 



the applied external field can be tracked down to the fact that one tries to model the effect 
of the battery within the calculation. 

The appropriate form of the applied field is available from considerations originally made 
in the linear response regime^. To have a system infinite along the direction of flow and 
characterized by a finite drop in external potential, AV, one needs to consider the large-time 
limit of a field 

E ext (x, t) = ^(0(x + ut) - 6(x - ut)), (2) 

where 9{x) is the unit step function, x is the direction of current flow and u is the desired 
speed of the front between the region with an applied field and that without field. This 
represents a situation that has a constant drop AV in the potential at all times t and, 
for sufficiently large u 3> vp, it produces a uniform current even in the case of interacting 
electrons, i.e. the field is not screened out. The steady state is obtained in the t — > oo limit 
where the drop AV persists, but the external field in any finite part of system is zero. (e.g. 
E(x = 0,t) = AV/(ut) — > 0). The initial localization of the field-containing region leaves 
no long-time signature other than the steady current that flows, and associated changes in 
the electronic structure such as density or induced potential. 

The above considerations of the steady state show that the Hamiltonian with zero applied 
field must be used when constructing the density operator within the maximum entropy 
ansatz. Similarly, we believe that the use of a finite external field together with the maximum 
entropy prescription is simply incorrect, and its application in other calculations should be 
reconsidered. The induced field will, however, appear in the calculation as a consequence of 
the current constraint applied to the density matrix. 

The second outcome of this observation concerns the steady-state character of the sys- 
tem. Once we accept that the Hamiltonian present is that without the applied field, the 
stationarity of the statistical density operator requires^ 



p,H 



(3) 



This identity has to be included when performing the constrained search for the operator p. 
At the same time, it should be clear that this condition can be fulfilled only for a system 
infinite along the direction of the current flow. (The only exception is a system of perfect 
translational invariance and finite periodic boundary conditions. As this represents a very 
special and highly non-generic case - an arbitrarily small perturbing potential spoils the 
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perfect translational invariance and therefore the ability of the system to carry current - we 
will not be concerned with it in our further discussion.) 

The condition (jHJ) has been shown to follow also from a time-evolution point of view by 
HershfielcU^. In this work it is also correctly pointed out that, as opposed to the equilibrium 
expectation value, non-equilibrium systems are characterized by an effective Hamiltonian 
that enters the statistical density operator, which is different from the true physical Hamil- 
tonian characterizing the time evolution or time-correlations in the system. This, as we 
will see, significantly complicates the formal development of the theory for interacting non- 
equilibrium steady-state systems. 



III. THE INVARIANT CURRENT APPROACH AND THE INDUCED DROP IN 
THE POTENTIAL 

In the invariant- current approach the constraint that keeps the system out of equilibrium 
is chosen to be the current at a particular point, xq 

I(x ) = JdS-Ti pj(x ,y,z) , (4) 

where j(x, y, z) is the operator of the current density at r = (x, y, z). As we have shown2&, 
this requirement, together with the steady-state restriction (jH)), leads to the statistical op- 
erator 

p = exp{n- P(H - fiN - AT )}. (5) 
where the operator 1° is the invariant part of the current operato r 30 ! 31 

i r T - 

lim — / I(t)dt, (6) 



T^oo 2T 



T 



which is independent of the position of current measurement xq. This arises from the fact 
that for a stationary density matrix the current fulfills the continuity equation V-j = —n = 0. 
The time- dependence of the operator I(t) is determined by the Hamiltonian of the system 
H which, similarly to the case of scattering field operators in Hershfield's worki£, hinders 
the development of the theory for interacting electrons. 

At the mean-field level of approximation, the one-particle density matrix can be used to 
obtain all required quantities. Instead of the many-body Hamiltonian we need to consider 
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n[x, x 



the effective one-particle Hamiltonian given bj»* 

h = -±% + V(x)-iJL-AI (x), (7) 

where V(x) is the self-consistent potential, \i chemical potential, A a Lagrange multiplier 
belonging to the current constraint and J° the invariant current operator. Even though 
we eventually obtain the one-particle density matrix that closely resembles the usual local 
Fermi distribution, its derivation is non-trivial. The problem arises from the fact that our 
system cannot, strictly speaking, be obtained as a limit of a finite one. First, the finite 
drop in the potential, A0, makes it impossible to consider periodic boundary conditions, 
and, second, the existence of a non-zero current flow hinders the construction of hard walls 
placed at finite, but large distances at right and left, as used by Adawi^S. Similarly, the 
use of periodic boundary conditions, as implicitly appear in some current-constraint-based 
treatments^^^, is not consistent with the possible existence of an overall drop in the 
electrostatic potential and the non-zero current. For these reasons we give its detailed 
derivation in the Appendix [XJ The resulting one-particle density matrix is 

+ 1 

where XE,a{ x ) are energy- normalized states that diagonalize the Hamiltonian (JJJ) and E a (E) 
are the corresponding eigenvalues. The latter can be expressed in the basis of scattering 
states of the physical Hamiltonian (i.e. without the term containing the current operator in 
Eq. ©) 

Xe,ol{x) = S ^<pE,v{ X ) U vA E )i ( 9 ) 

V 

where 0£,„(x), v — R,L represent right- or left-going scattering states at the energy E given 
asymptotically as 

» 0) = "y= te ^> ( 10 ) 
<PeA* « 0) = -^=te- ikx , (11) 

where k = y/2E and q = \/2(E + A0) are the wave-vectors on the far left and far right 
respectively and t and t are transmission amplitudes for right-going and left-going electrons 
respectively, (t and t depend on the energy but we will not write this dependence explicitly.) 
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^ _ V- [ , f XeAx)XeA x ') rs x 



We note that the states XE,a can be labeled with the energy because the effective Hamiltonian 
commutes with the physical one, as required by the stationarity condition Q. From this 
and a glance at Eq. (j7J) it follows that the states XE,a are simultaneously eigenstates of the 
invariant current operator 7°. Therefore the second index differentiates between energy- 
degenerate states which, in the simplest case of a single ID channel considered here, attains 
two different values, '+' and depending on the sign of the invariant current eigenvalue of 
the respective state. In Appendix [B] we show that the states Xe,+ and Xe.- transform one 
into another under the time-reversal T, 



Xe, 



e^XE, 



(12) 



The latter relation has an important consequence for the induced change in the density in 
the linear regime, as we will discuss below. 

We have already indicated that the scattering-states representation plays a fundamental 
role not only in the LRS approach but also in the current-constraint schemes. It is therefore 
useful to express the current operator in the scattering states representation^ 



d^jx « 0) f 



2tt 



dE 

dl^ v (x » 0) 
dE 



1 - \r\ 2 —r*i 
-rt* 



±\t\ 2 



27T 



dl v ^(x « 0) 
dE : 



(13) 
(14) 



where S V)V is the scattering matrix 



Sun 



r t 
t f 



(15) 



and the last equality in Eq. (|14j) follows from the unitarity of S u ^ (current conservation) 
and t, r and t, r are the transmission and reflection amplitudes for right- and left- going 
electrons. Similarly to the latter, we do not write explicitly the energy dependence of the 
S'-matrix or the current operator matrix. Since the invariant current operator 7° in (|7J) is 
related to the matrix of the current operator multiplied by a delta-function of energy (see 
Appendix IXj) . 



I°(E,a;E',p) x A 



dl 



a,/3 



dE 



5(E - E') x A, 



(16) 



the states XE,a{x) automatically diagonalize the kinetic and potential energy terms in 
Eq. ((7|). To diagonalize the complete effective Hamiltonian they need to diagonalize also 



the current operator which is given by Eq. The unitary transform that achieves this has 
been found before^ 



r*t/\t\ r*t/\t\ 



v/2(Hti) v^+i*i) 
1*1—1 iti+i 



_ y/2{l-\t\) y/2(l+\t\) 

with the corresponding effective eigenvalues 

E a (E) = E±A\t\. 



(17) 



(18) 



The calculation of the current is then straightforward using the one-particle density matrix, 



(/(*)> = E j dE - 



ME a (E)-p) 



t dI U)V (x) 



(19) 

7 + 1 

Using this expression we find the dependence of the current on the renormalized Lagrange 
multiplier A. 

For comparison with other approaches as well as experiments, we also need the depen- 
dence of the current on the drop in the electrostatic potential, A0, which can be obtained 
from the induced change in the density Sn(x) via the expression 



(20) 



4>{x) = I v(x,x')Sn(x')dx' 
A(f) 



-oo — 



-oo 



where v(x,x') is an appropriate effective electron-electron interaction. Alternatively, A<p 
can be obtained from the local neutrality conditions in the right and left electrodes^ for the 
self-consistently determined scattering matrix (hence the link to the resistivity dipoles based 
formula). The local neutrality conditions demand 



n(x « 0) + n B = 
-n(x » 0) + n B = 0, 



(21) 
(22) 



where n B is the positive charge density of the background, assumed to be the same in both 
electrodes for simplicity. Subtracting these two we get a single condition 



An = n(x « 0) - n(x » 0) = 



(23) 



This identity alone can be used to determine the drop in the potential A0, for given values 
of A and /i. The density in the asymptotic regions is obtained in a way similar to the total 
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current: we first find the expressions for the matrix of the density in the asymptotic regions 
within the scattering states representation 

dn vr] {x « 0) Ire 1 ct c 1 x x 1 dl uv {x « 0) 

dE = fcrife^' 1 + = - * dE (24) 

dn^(a; >> 0) 1 1 t _ 1 l di^Qr >> 0) , . 
dE - 2^ 2 ^' 2 + 2^> 2 ^ " ^ A ' 2 + g dE • (25) 



Second, we express the local neutrality condition (|23|) using the density matrix 

dn v ^(x « 0) dn v ^{x » 0) 



y / ^ 1 — c/t 

^ / e f)(E a (E)-p) + 1 



£V = 0. (26) 



a, u,r] 

To obtain algebraic results we need to restrict our treatment to the j3 — > oo limit (corre- 
sponding to the density matrix with the minimal energy for given constraints). Under these 
circumstances the effective Fermi distribution takes the form 

R(p rF | x « ^O* - ^ + + *«^0* - ^ - 4*(*0I). (27) 

where the two step functions #() correspond to the positive and negative eigenvalues in 
Eq. ED Assuming A > 0, we see that starting from E > /i 2 = \i — v4|t(/i 2 )| and up to 
E < fii = fi + A\t(ni)\, only one of the two degenerate states XE,a will be occupied^. 

For E < fj,2 the contributions to the density clearly can not depend on the unitary rotation 
between the scattering states and we therefore have 

£ *^{z « 0) = £ « = ± (28) 

V %^(x»0) = — . (29) 
^ dE V ' irq K ' 

a 

The contribution from the singly occupied state XE,a=i{%) we find from Eqs. ()24I25I17|) with 
a surprisingly simple result 

^«°) = S^(*«»)^ = SE. (») 

^»0) = -L. (31) 
dE K J 2nq K J 

Finally combining Eqs. (|27I29I28I EUJ EH) in the Eq. fl2EJ we get 

^ 2 dE dE f 112 dE r 1 dE 

o 2nk 7„ A ^ vrg J ^ 2nq 
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The last equation can be satisfied only for A0 = 0. This means that the invariant-current 
scheme with the off-diagonal elements at the same energy retained leads to no induced 
potential drop and therefore its applicability to common nano-contacts is doubtful, even in 
a const ant- current experiment. 

Qualitatively, this result is a consequence of the fact that the contribution to the density 
per-energy of each doubly degenerate state XE,a, dn atCt /dE, is almost the same^ far right 
and far left (see Eqs. ([3U|3ip ). If we assume that the right-electrode, having the same 
background charge density as the left one, has the bottom of its local density of states at 
E = — A0, below that of the left electrode {E = 0), occupying this portion of the energy 
spectrum (E < 0) will partially compensate for the background charge in the right electrode. 
Once the states below the bottom of the left electrode are filled, adding each electron into 
next state XE,a will contribute in both electrodes almost equally so that when attempting 
to compensate for the background charge of the left electrode we will inevitably overload 
the right electrode or, when neutralizing the right electrode there will not yet be enough 
electronic charge in the left one. (We would like to stress that the above reached conclusions 
are valid for (3 — > oo limit, finite A and a regime, in which the equations for /ii and fi2 have 
a unique solution. The latter fails to be fulfilled if either fii or /i 2 approaches a resonant 
energy level of the potential V(x).) 

An even more surprising result appears in the linear- response regime, i.e. when AI° << fi. 
Under these circumstances the induced change in the density 5n(x), 

Sn(x) = n(x, x) \ A - n(x, x)\ A=0 , 

where n(x, x)\ A is the diagonal of the density matrix given in the Eq. (JHJ) for a given (small) 
value of the renormalized Lagrange multiplier A, is simply given by 

5n(x) = 2A\t(p)\ (|x+(*)| 2 - |x^)| 2 ) + 0(\A\ 2 ). (33) 

Using the time-reversal character of the states XE,a, Eq. (jl2)l . we immediately obtain that 
5n(x) = 0(\A\ 2 ), i.e. there is no change of the density in the linear regime. Using this result 
in the formula for the induced drop, Eq. (J2"U|) . we once again confirm the above obtained 
result (valid even for finite A) of zero induced drop in electro-static potential. 
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IV. THE INVARIANT CURRENT APPROACH WITH DECOHERENCE 



The situation is quite different when the off-diagonals of the invariant current operator 
are dropped. This makes the density matrix (jSJ) identical to that given by Ngi^. The 
elimination of the off-diagonals can be made physically plausible from the assumption of 
phase independence between the right- and left- going scattering states. Namely if we 
ascribe an independent fluctuating phase e tdv to each of these two, the off-diagonal matrix 
elements of the one-particle density matrix, Eq. JSJ), are easily seen to be averaged to zero. 
The effect of the coupling to the environment is to suppress, through averaging over a 
fluctuating phase, any mixing of left- and right- going states in the pure states that are 
summed in the density matrix. In effect, these pure states are constrained to be either left- 
or right- going. This means that the off-diagonal elements of the current operator no longer 
play any role in the determination of the steady-state density matrix. 

To recast the above considerations into a more formal language we need to include a 
description of the effect of the environment in the Green's-function-based derivation of the 
one-particle density matrix given in Appendix (jXJ. The Green's function, that directly leads 
to the one-particle density matrix, is defined in terms of the time-ordered product of the 
field operators, 

G(2,l) = -Tr[pT{^(2)^(i)}], (34) 

where p is the many-body density matrix that specifies occupations of the states of the 
whole system, i.e. the degrees of freedom of the environment that cause the above discussed 
fluctuating phases. Rewriting the Green's function in terms of creation and annihilation 
operators of electrons in the scattering states, 

G(2, 1) = - J dEdE' ^,,(x)^,,(x')Tr [pT{c^(^4, , (35) 

i/,rj 

we observe that the Green's function will have only the diagonal elements with respect to 
the v,T] = R/L index, i.e. 



Tr 



/df) df) 
^£^p e e ie »e-^T{c E At2)J E , v m 



Su, v , (36) 



where p e is a many-electron density operator already without the environment's degrees of 
freedom which were effectively taken into account via averaging over the phases of the right- 
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and left- going scattering states. From this it follows that the equation of motion for this 
Green's function, Eq. flA8|) . is already in a diagonal form and so is the one-particle density 
matrix n(x,x f ), Eq. (jA17l) . 

Anticipating the sources of decoherence and therefore the irrelevance of the invariant cur- 
rent operator's off-diagonal elements, we can proceed rather straightforwardly in a derivation 
of Ng's^ as well as our former results^ for the theory with decoupled right- and left- going 
states. Since the current matrix (|T3*j) is effectively in a diagonal form, we can obtain the re- 
sults considering the transform U U)OC = 5 U)a , so that the states XE,a(x), entering the one-body 
density matrix (jHJ), are directly the scattering states. The effective eigenvalues prescribing 
the occupations are then given by 

E a {E) = E±A\t{E)\ 2 . (37) 

What comes as an essential difference, as compared to the case when the decoherence is not 
accounted for, is that the contribution to the density from the singly-occupied states now 
gives 

as can be easily seen by inspection of the expression (|24|) . (J25|) . (|13|). Clearly, the contribu- 
tions to the local densities far right and far left are now significantly different and therefore 
the qualitative argument for zero induced drop in the potential does not apply. Using expres- 
sions (|3~% j) .(|3T H) in the local neutrality condition (j2HJ) eventually leads to the linear response 
result 

A0 = 2i|t(/i)| 2 |r(/i)| 2 . (40) 
Using equations (fT§|) . (|2"7j) and (fT3*|) we obtain the current in the linear response 

(I) = 2A\t(^ l -^f. (41) 
and therefore combining the last two expressions we arrive at the well known result 

±_ = UML (42 ) 

for the so called 4-point conductance, a result already mentioned also within current con- 
straining scheme o 15 i 26 and in agreement with the LRS and therefore it is also compatible 
with the linear regime of the Landauer-Biittiker formalism. 
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To summarize this section, incorporating a decoherence caused by many-body interac- 
tions, the invariant current approach does lead to a non-zero induced drop and the resulting 
conductance is identical to that obtained within NEGF formalism as well as the LRS. We 
should stress that the way we have incorporated the decoherence contains assumptions re- 
garding the coupling of the states to the environment and therefore can not be regarded as a 
truly ab initio approach. The second important outcome is that formulating the maximum- 
entropy current-constrained scheme that starts at the non-interacting or mean-field self- 
consistent field level will not give physically meaningful results. 



V. THE UNIFORM-CURRENT THEORY 

In the uniform-current theor y 2 ?' 24 ' 2 ^ one takes the constraint to make the current uniform 
throughout the system 

/ = /(*) = (?(*)>• (43) 

This is achieved by introducing a continuous Lagrange multiplier A(x) and the corresponding 
DM takes the form 

p uc = exp ^Q uc — (3 — jiN — j dxA(x)I(x)\ j (44) 

The function A(x) has to be found such that the constraint f)43j) holds. The essential 
difference from the invariant-current scheme is the fact that this ansatz results in a density 
matrix which evidently does not commute with the physical Hamiltonian 



P UC ,H 



? 0. (45) 



As a consequence, even though the current is uniform at some instant, (and therefore from 
the continuity equation the density is momentarily stationary), it will in general change at 
later times together with many other averages computed using the p uc . 

Apart from this objection, one can show that the scheme is equivalent to a equilibrium- 
like calculation with some effective potential, using a gauge transformation. Similarly to 
the invariant current theory, it emerges that the induced drop in electrostatic potential is 
zero. We will demonstrate this at the mean-field level of approximation only, even though 
it presents no complication in this case to prove it for fully interacting electrons as well. 
We firstly reformulate the many-particle problem of the full density matrix p uc into that 
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of one-body density matrix, in a way completely analogous to the previous section and the 
formalism given in the Appendix [X] Eventually we will be concerned with the mean-field 
single-particle Hamiltonian of the form 

h = -\d 2 x + V(x) + f d 3 x'A(x')i(x') = \ (-id x + A(x)f + V{x) - \A\x), (46) 
z j z z 

where A(r) is the same Lagrange multiplier as used in the many-particle density matrix, the 
current operator is 

W) = -\ (dJ(x - x') + 6{x - x')d x ) , (47) 

where the subscript x reminds that the operator operates on functions of x. In this way 
the non-equilibrium problem has been transformed into a complex (so that current flow is 
possible in an effectively equilibrium-like system) but Hermitian eigenvalue problem on the 
whole space x G (—00, 00) (so that the eigenstates form a continuum) 

hxE,a( x ) = ExeA x )- ( 48 ) 

The final one-particle density matrix is given by analogy with the previous section by the 
expression 

n(x, x ) = ^ / dE e p {E _ ll) + 1 ■ (49) 

a 

However, the complex character of the problem can be removed by a simple gauge transfor- 
mation 

XEA x )=e- irA{x ' )dx '<pEA x )i (50) 
where the transformed states 4>E,a are eigenstates of a real Hamiltonian 

h' = - l -dl + V{x)- l -A\x). (51) 

From this follows that the one-body density matrix 

n(x,x') = W d Ee^rAO^tE0^^1 e ^r'Am^ (52) 

a 

gives non-zero current only through the gauge-factors. The current can be therefore easily 
evaluated to give 

(I(xo)) — / dxdx'5(x — x')I x (xo)n(x, x') = A(xo)n(xo), (53) 
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where n{x) is the electronic density. We can now fulfill the requirement of uniformity of the 
current by giving the Lagrange multiplier as 

A(x) = -L-, (54) 
n[x) 

which makes the similar results of the previous work by Kosov^ completely general (the 
same result holds for interacting electrons since the gauge transform argument does not 
depend on the interactions). Using the result (|54jl within the effective Hamiltonian (|51j) 
gives a simple closed set of equations to be solved. 

Finally we turn to the analysis of the induced drop in the electrostatic potential which, 
similarly to the case of the invariant current theory, is identically zero. To show this we note 
that the self-consistent potential V(x) used to determine the wavefunctions and therefore via 
Eq. lJH2]l the density does not have any finite drop as it corresponds to a fictitious equilibrium 
system (the effective Fermi function in (JH2|) depends only on the energy of the single-particle 
wavefunction <pE,a{x)-) for which the drop must be clearly zero. This observation is not 
affected by the presence of the last term in the effective Hamiltonian (j51|) since in the case 
of identical electrodes 

A(x oo) - A(x -> -oo) = / (l/n B - l/n B ) = 0. (55) 
and the contribution so the drop in potential is zero. 



VI. CONCLUSIONS 



In conclusion, we have shown how the maximum-entropy formalism can be applied for 
non-equilibrium steady states within the framework of the single-particle approximation. We 
have presented three different approaches within a common formalism: (1) the invariant- 
current constraint, (2) the invariant-current constraint without the off-diagonals and (3) the 
uniform-current constraint scheme. For these we have rigorously derived the one-particle 
density matrix for an infinite system that cannot be realized as a limit of a finite system. 
Subsequently we have obtained the expression for the electrostatic drop between electrodes 
within the linear response regime. We have shown that (1) and (3) give zero induced drop 
in electrostatic potential which is not compatible with actual current- carrying situations in 
nano-j unctions. In the case (2) we have shown that by removing the off-diagonal elements 
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of the current operator, the induced drop is nonzero and in fact, the results in the linear 
response are identical to the LRS. We can view the three different schemes analyzed in the 
two last sections from a more general perspective. They represent a variational prescription 
for the search of the non-equilibrium state with a given average of total current. The key 
difference between them is the part of current operator that is being used for the constraint. 
The uniform-current theory takes the whole, unmodified current operator; the invariant- 
current scheme (through the stationarity requirement) removes the off-diagonals between 
states of different energy; and finally there is the form of invariant- current theory in which 
the off-diagonal elements within the scattering states basis set representation are removed 
by decoherence. Interestingly, only the last gives a electronic density which results in a 
nonzero induced drop in the potential. This observation rises a question whether there is 
something inherently wrong with demanding the current to be fixed. The theory with no 
off-diagonal elements of the current operator shows that the explanation for the problems 
with the current-constraining schemes arise from the fact that many-particle interactions 
or a certain source of decoherence is essential for the density matrix to give a physically 
meaningful results. 
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APPENDIX A: DERIVATION OF THE ONE-PARTICLE DENSITY MATRIX 

To derive the form of the one-particle density matrix we first define an auxiliary Green's 
function, in a close analogy with the Matsubara technique^, 

G(2,l) = -Tr[pT{>(2),^(l)}], (Al) 
^t(!) = e I<T1 ^(x 1 )e- kT \ij(2) = e &T ^(x 2 )e- kr2 (A2) 
K = H-/2N-AI , (A3) 
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where ift(x) and ip(x) are the field operators for electrons. The fictitious time dynamics is 
given by the effective many-particle Hamiltonian K, such that one can employ the similarity 
between the unitary time evolution operator and the density matrix p The Green's 
function defined in this way satisfy the equation of motion 

\-\% + V ^ AT°(z) + d T ^j G(x, x') = -5(x - x'), (A4) 

with the fermionic boundary condition 

G(t = 0,r') = -G(t = /3,t'). (A5) 



We note that while the many-body effective Hamiltonian K in (|A3|) is infinite and therefore 
just a formal expression, the one-body Hamiltonian present in the equation of motion for the 
Green's function (IA4J) is finite, for it represents energy per particle. The scattering states 
4>e,v(x) (jTDjl. (fTTj) diagonalize the Hamiltonian — \d"l + V(x) and therefore leave the elements 
of 1° in a block-diagonal (diagonal with respect to the energy) form given by2& 

I°, V {E,E') x A = lim -j- / dt(0 B)i/ |/>,t)|0 E! A x A (A6) 

= « 0) x 5{E - E') x A, (A7) 

dE 

where we have used the renormalization of the current-related thermodynamic parameter A 
into A = limr^oo irA/T. The matrix dI v<n jdE is the matrix of current operator at energy E 
given by Eq. ([13]) . The equation of motion for G(r) in the scattering-states representation 
now takes the form 



(E - p)5 vy - A^f + d T 



Ge'Ut, t') = -S(t - r')5(E - E")S^, (A8) 



E 

u' 

where 

GTA^r') = («^|G(x,:zOI<W- (A9) 



It is now clear that to solve ()A8|) we finally need to diagonalize the current matrix, i.e. 
employ the unitary transformation U a>v discussed in the section (jlll|) . The equation of 
motion after this transform takes a simple from 

(E a - d T ) G a E %(r) = -5(t - T')5 a ,p5(E - E'). (A10) 
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We will drop the energy E, E' and state a, (3 indices for a moment since all the following 
manipulations are diagonal with respect to these. Regarding the time, the Green's function 
must obey the fermionic boundary condition 

G(0) = -G{-/3). (All) 

The general solution of Eq. (jAlOj) is clearly 

f Ce-^ T+ V -(3<r<r' 
G = I . (A12) 

[ -Ce~ Er t 1 < r < 

Integrating the equation of motion over f T ,^$ we see that G(r) must have a unit step 
discontinuity at r = r' so that we can fix the constant C 

Ce -E(T> +I 3) _ (_c e -^') = 1 (A13) 

and therefore 



C = — (A14) 



The density matrix is then given by 



Et> e -E(T+f3) f 3<T<T > 

G(T,T') = -—rx\ P (A15) 

1 + e~ E - e ~ ET t' < r < 



and the density matrix through the prescription is 



-EQ- 

= G( r = IT, 0) = -g— ■ (A16) 



Restoring all the indices we finally have 

,XEA X )X*E,a( X ') 



n(x, x') 



G(x, r = 0-; x f , t' = 0) = Y f dg ^W) < (A17) 



Which is the result (JBJ) given in the section (IIII|) . 

The above given derivation also shows that our approach does not depend on the choice 
of the normalization. If we had chosen the scattering states normalized to a delta-function 
of k instead of E, the current matrix elements (|13|) would have to be multiplied by k. 
However, the delta-function of energy in ()13|) is unchanged, since it comes from the general 
considerations of stationarity of the ensemble^, so that we have kS(E — E') = 8(k — k'). 
Using this in the equation of motion for G we see, that the final result is the same, i.e. it 
is again the matrix elements of current in the energy-normalized states that appear in the 
effective dispersion relation. 
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APPENDIX B: TIME-REVERSAL SYMMETRY AND THE MAXIMUM EN- 
TROPY STATES 



In this Appendix we show that the states x+( x ) an d X-( x ) are related by the time- reversal 
symmetry 

T X+ (x)=xU x ) = ei *X-(x). (Bl) 

Let h° = h + AI° be the physical one-particle Hamiltonian in the Eq. (jJJ). Next, let 
are the single-particle states that diagonalize both h° and 1°. We will use T for the 
time reversal operator which is simply a complex conjugation. We have 

h° \ X+ ) = e \ X+ ) (B2) 

h°T \ X+ ) = eT \ X+ ) , since [H°, T] = (B3) 

1° \x+) = +^ \X+) (B4) 

I°T \ X+ ) = -iT | X+ ) , since TI° = -I°T (B5) 

The last property of the current operator is true for representation of operators and wave- 
functions in a real space, where complex conjugation of the current operator changes it's 
sign. Identical statements hold for the left- current- carrying maximum-entropy state |%_) 
only with a reversed sign of the current eigenvalue i. 

We know that \x+) an d \x~) are degenerate eigenstates with an eigenvalue e. From 
Eq. (jB3|) we see that T\x+) is also a state degenerate with them, and from Eq. (jB5|) that it 
has the same current eigenvalue as |%_). Since this exhausts the possible degeneracy (2-fold 
in ID), the only possibility is that 

T\ X +) = e«\x-), 

where is an arbitrary phase factor. This therefore shows that \x+) an d \x~) are related 
by time-reversal symmetry. 
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